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Summary 

• Photoperiodic flowering is a major factor determining crop performance and is controlled 
by interactions between environmental signals and the circadian clock. We proposed Hvluxl , 
an ortholog of the Arabidopsis circadian gene LUX ARRHYTHMO, as a candidate underlying 
the early maturity 10 (eamW) locus in barley (Hordeum vulgare L). 

• The link between eam10 and Hvluxl was discovered using high-throughput sequencing of 
enriched libraries and segregation analysis. We conducted functional, phylogenetic, and diver- 
sity studies of eam10 and HvLUXI to understand the genetic control of photoperiod response 
in barley and to characterize the evolution of LUX-Wke genes within barley and across mono- 
cots and eudicots. 

• We demonstrate that eam10 causes circadian defects and interacts with the photoperiod 
response gene Ppd-H1 to accelerate flowering under long and short days. The results of 
phylogenetic and diversity analyses indicate that HvLUXI was under purifying selection, 
duplicated at the base of the grass clade, and diverged independently of LUX-Wke genes in 
other plant lineages. 

• Taken together, these findings contribute to improved understanding of the barley circa- 
dian clock, its interaction with the photoperiod pathway, and evolution of circadian systems in 
barley and across monocots and eudicots. 
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Introduction 

Barley {Hordeum vulgare L.) is an important crop cultivated in a 
wide range of environments and has recently emerged as a model 
species for Triticeae as a result of its diploid nature and abun- 
dance of genomic resources (International Barley Sequencing 
Consortium (IBSC), 2012). Genetic variation in photoperiod 
response was crucial for the successful expansion of barley cultiva- 
tion from its origin in the Fertile Crescent to northern latitudes 
(Turner etal, 2005; Jones etai, 2008). Barley is a facultative 
long-day plant; long days (LDs) promote flowering in spring, 
while short days (SDs) delay reproductive development. Flower- 
ing under LDs is controlled by the major photoperiod response 
PHOTOPERIOD 1 gene {Ppd-Hl) (Turner etal., 2005). 
Ppd-Hl is homologous to PSEUDO-RESPONSE REGULATOR 
(PRR) genes implicated in the circadian clock of the model spe- 
cies Arabidopsis thaliana (hereafter Arabidopsis, Turner etal., 
2005). The dominant allele of Ppd-Hl is prevalent in the wild 
barley progenitor H. vulgare spp. spontaneum and in Mediterra- 
nean cultivated barley genotypes and accelerates flowering under 
LDs, as an adaptation to short growing seasons. A natural muta- 
tion in the conserved CCT domain of Ppd-Hl causes a reduced 
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response to LDs and was selected for adaptation to long growing 
seasons (Turner etal, 2005; von Korff etal, 2006, 2010; Jones 
etal, 2008; Wang etal, 2010). Under LD conditions, Ppd-Hl 
up-regulates HvFTl, which is the barley counterpart of the Ara- 
bidopsis 'florigen' FLOWERING LOCUS T(ET). Delayed flow- 
ering in genotypes with a mutated ppd-Hl allele was associated 
with reduced expression of HvFTl (Turner etal, 2005; Campoli 
etal, 2012a). 

Barley early maturity mutants {earn) with a reduced or no 
response to photoperiod have been isolated and used in Swedish, 
Australian, and South American breeding programs to adapt cul- 
tivars to short growing seasons (Laurie etal, 1995; Lundqvist, 
2009; Zakhrabekova etal, 2012). These earn mutants have been 
introgressed into the spring barley cv Bowman (Druka etal, 
2011) and represent a valuable resource to decipher the genetic 
control of photoperiod response in the model crop barley. 
Recently, HvELF3, a homolog of the Arabidopsis circadian clock 
regulator EARLY FLOWERLNG 3 (ELF3), was identified as a 
gene underlying the eam8 quantitative trait locus (QTL), which 
causes a day- neutral early flowering phenotype (Faure etal, 
2012; Zakhrabekova etal, 2012). Interestingly, the mutation in 
HvELF3 caused an up-regulation of Ppd-Hl and the downstream 
HvFTl under noninductive SD conditions (Faure etal, 2012). 
In addition, the Hvelf3 mutants were severely compromised in 
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the expression of clock oscillator and output genes. This study 
suggested that circadian clock homologs play an important role 
in the control of flowering in barley. Circadian clocks synchro- 
nize biological processes with the diurnal cycle, using molecular 
mechanisms that include interlocked transcriptional feedback 
loops. In Arabidopsis, the circadian clock is composed of three 
negative feedback loops: the inhibition of evening complex (EC) 
genes ELF3, EARLY FLOWERING 4 (ELF4), and LUX 
ARRHYTHMO {LUX, also known as PHYTOCLOCK1) by the 
rise of CIRCADIAN CLOCK ASSOCIATED1 (CCA1) and 
LATE ELONGATED HYPOCOTYL (LHY) late at night; the 
inhibition of PRR genes by the EC early at night; and the inhibi- 
tion of LHY/CCA1 by TIMING OF CAB EXPRESSION1 
(TOC1) in the morning (Kolmos etal., 2009; Huang etai, 
2012; Nagel & Kay, 2012; Pokhilko etal., 2012). Furthermore, 
the evening-expressed GIGANTEA (GI) protein was modeled as 
a negative regulator of the EC, which in turn inhibits TOC1 
expression (Herrero etal., 2012; Pokhilko etal., 2012). 

Campoli etal. (2012b) have shown that circadian clock genes 
are structurally conserved between barley and Arabidopsis, and 
their circadian expression patterns suggested conserved functions. 
However, phylogenetic analyses revealed that duplications/dele- 
tions of clock genes occurred throughout the evolution of eudi- 
cots and monocots. For example, the ancient three PRR clades 
expanded independently in monocots and eudicots, supposedly 
as a result of paleoduplications (Takata etal., 2010; Campoli 
etal, 2012b). In monocots, these events gave rise to two paralo- 
gous pairs of PRR genes, termed PRR37I73 and PRR59/95. In 
this context, it is noteworthy that /Wy7orthologs in monocots, 
PPD1 in barley and wheat (Turner etal, 2005; Beales etal, 
2007) and SbPRR37 in sorghum {Sorghum bicolor) (Murphy 
etal, 2011), are major determinants of photoperiod sensitivity 
and flowering time, while natural variation in PRR genes in 
Arabidopsis did not have a strong effect on flowering time 
(Ehrenreich etal, 2009). These results indicated that the molecu- 
lar function of genes participating in the circadian clock and pho- 
toperiod response diverged and specialized in a lineage-specific 
manner. Because of the prominent role of photoperiod insensitiv- 
ity in breeding, it is essential to identify the genetic components 
controlling this pathway in temperate cereals. 

We present the characterization of the barley eamlO mutation, 
which causes an early-flowering phenotype under both SDs and 
LDs. eamlO was described as an X-ray-induced mutation in 
Super Precoz (2H) of unknown parentage (Favret & Ryan, 1966; 
Gallagher etal, 1991; Gallagher & Franckowiak, 1997) and was 
mapped onto chromosome 3HL (Borner etal, 2002). Using a 
high-throughput sequencing approach, we identified the candi- 
date gene underlying the earn 10 locus in barley as HvLUXl, a 
barley ortholog of the Arabidopsis circadian clock regulator LUX. 
We demonstrated that eamlO caused circadian defects and inter- 
acted with the major barley photoperiod response gene Ppd-Hl 
to accelerate flowering under LD and SD conditions. Based on 
the results of phylogenetic and diversity analyses, we conclude 
that HvLUXl was under purifying selection, duplicated at the 
base of the grass clade, and diverged independently of the L UX- 
like genes in other plant lineages. 
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Materials and Methods 

Plant material and growth conditions 

Flowering time (days to emergence of the main spike awns) of 
the barley spring cv Bowman and three Bowman-derived intro- 
gression lines Bowman(/))(/-ffl), Bowman{eamlO), and Bow- 
mzn{Ppd-Hl + eamlO) (kindly provided by R. Waugh, James 
Hutton Institute, and by David Laurie, John Innes Centre, UK) 
was recorded for 15-18 plants per genotype. To score flowering, 
plants were grown in soil in a glasshouse under SDs (8 h light, 
20°C: 16 h dark, 18°C) and LDs (16 h light, 20°C : 8 h dark, 
18°C). To investigate expression levels of HvFTl, plants were 
sampled 2h before light-off after 14 and 28 d under LDs and 
SDs, respectively. Diurnal and circadian expression of core clock 
and flowering time genes was tested under SD and free-running 
conditions in Bowman and Bowman(^w70). Plants were grown 
in soil in a controlled-environment growth chamber (CEGC). 
After 21 d under SDs, leaf material was harvested every 2 h for a 
total of 24 h from the start of the light period (time point TO). 
Night samples (T10 to T22) were collected in the dark. After 
SDs, plants were released into continuous light and constant tem- 
perature (20°C; LL) and sampled every 2 h for 48 h, starting after 
8 h of continuous light (T8). Two biological replicates, compris- 
ing the second youngest leaves of three independent plants, were 
analyzed. All samples were immediately frozen in liquid nitrogen 
and stored at — 80°C until processed. 

Meristem development was scored in a separate experiment 
with Bowman and ^>ovfmm{eaml0) grown in soil in a CEGC. 
The main stem of three plants per genotype was dissected starting 
7 d after germination every 2-3 d under LDs and every 3-4 d 
under SDs until flowering. After day 79 under SDs the experi- 
ment was stopped when the shoot apical meristems of the 
remaining plants had died. Meristem development was scored 
following the Waddington scale (Waddington etal, 1983). To 
generate developmental series for the analysis of HvFTl expres- 
sion in the leaf and to link it to specific meristem stages, two bio- 
logical replicates were harvested 2 h before light-off at every 
developmental stage under both photoperiods. 

To investigate natural diversity of a candidate gene, a set of 88 
wild {H. vulgare ssp. spontaneum and H. agriocrithon) and culti- 
vated barley genotypes was selected from germplasm collections 
of the Max Planck Institute of Plant Breeding Research in 
Cologne, Germany (Badr etal, 2000) and the Barley IK collec- 
tion (Hiibner etal, 2009; Supporting Information, Table S2). 

Gene expression analysis 

Gene expression was analyzed in leaf samples harvested from 
Bowman, J5owman{Ppd-Hl), Bowman(ftzm70), and Bowman 
{Ppd-Hl + earn 10) grown under LDs and SDs, from a develop- 
mental series of leaf samples from Bowman and Bowman 
{eamlO), and from diurnal and circadian sampling of Bowman 
and Bowman(ft2?«70) leaves. Total RNA extraction, cDNA syn- 
thesis, and quantitative reverse transcription polymerase chain 
reactions (qRT-PCRs) using gene-specific primers were 
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performed as explained in Campoli etal. (2012b). Expression of 
HvLUXl was analyzed using the following primer combinations: 
LUX_1077F 5 '-AATTCAGTCCACGGATGCTC-3' and LUX 
_1289R5'-CTTCACTTCAGCTCCCCTTG-3'. 

Identification and characterization of a candidate gene 
underlying eamW 

To identify a candidate gene for the earn 10 QTL, we sequenced a 
set of flowering-related genes from the genotypes Bowman and 
Howman(eamlO) using high-throughput sequencing of multi- 
plexed enriched libraries. A sample of pooled barcoded TruSeq 
libraries was prepared following the standard Illumina protocol 
with modifications (Methods SI). The pooled library was 
enriched for flowering-related genes using the SureSelect target 
enrichment system (Agilent Technologies, Boblingen, Germany) 
according to the manufacturer's recommendations with the fol- 
lowing modification. To alleviate biased PCR amplification of 
GC-rich regions, the post-enrichment amplification of the library 
was performed using the 'long denaturation' Phusion (Thermo 
Scientific, St Leon-Rot, Germany) protocol (Aird etal., 2011). 
The enriched library was paired-end-sequenced on a HiSeq 2000 
machine (Illumina, San Diego, CA, USA). 

The genes for targeted enrichment were selected by following 
two approaches. First, known flowering-related barley genes and 
gene families were directly extracted from NCBI GenBank 
(Table S3a). Second, additional barley genes were selected based 
on their homology with the genes of interest from 
Brachypodium distachyon as determined by the Blastx search in 
HarvEST assembly 35 (http://harvest.ucr.edu) (Table S3b). The 
nomenclature and annotation of Brachypodium genes were as 
reported by Higgins etal. (2010). Selected barley genes were used 
as a template to design a library of baits for the solution-based 
target enrichment approach and as a reference target genome for 
mapping in subsequent bioinformatics analyses. 

Read datasets were preprocessed using a set of filters (adapter 
trimming, quality and length, sequencing artefacts, contamina- 
tion) using the Galaxy server at Wageningen University, the 
Netherlands (http://galaxy.wur.nl; Methods S2) and mapped to 
the reference target genome as a single-end dataset using BWA 
0.5.9 with customized parameters: '-n' 0.03, '-o' 2, '-e' 10 (Li & 
Durbin, 2009a). Reads mapped to several locations and PCR 
duplicates were removed using 'view' and 'rmdup' functions of 
SAMtools 0.1.18 (Li etal, 2009b). Mean depth of coverage was 
estimated using BAMstats 1.25 (http://bamstats.sourceforge.net). 
Polymorphisms were discovered using 'IndelRealigner' and 'Uni- 
fiedGenotyper' algorithms of the GATK 2.1.3 package (DePristo 
etal, 2011). A set of stringent 'hard' filters was applied to the 
raw single nucleotide polymorphism (SNP) set to select reliable 
SNPs (Methods S3). 

Genes polymorphic between Bowman and Bowman(ftzw70) 
were located on barley linkage groups using GenomeZipper 
(Mayer etal, 201 1). The location of the eamlO QTL on the bar- 
ley consensus genetic map (OPA123-2008 consensus; Close 
etal, 2009) was extracted from the GrainGenesCMap browser 
(http://wheat.pw.usda.gov/cgi-bin/cmap). The position of the 
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MYB domain on HvLUXl was determined using InterProScan 
(Quevillon etal, 2005). 

Segregation analysis 

A set of 1002 BC1T2 plants derived from backcrossing Bowman 
(earn 10) to Bowman were grown in the field in Cologne in 2012 
(spring sowing). A total of 215 (21%) plants were scored as early 
flowering and seeds were harvested. Early-flowering BCiF 2: 3 lines 
were sown in the glasshouse under LDs. Leaf samples were har- 
vested and flowering time was scored. Barley genomic DNA was 
extracted using the BioSprint 96 kit (Qiagen, Hilden, Germany) 
according to the manufacturer's recommendations and quantified 
using Quant-iT™ PicoGreen assay (Invitrogen) measured using 
the Synergy™ 4 microplate reader (Biotek, Bad Friedrichshall, 
Germany). A full-length coding sequence (CDS) of HvLUXl was 
amplified using specific primers LUX_135F_T3 5'-aattaaccctcac 
taaagggTGGCGAGGGTAAGTTGATTC-3' LUX_1 108R_T7 
5'-gtaatacgactcactatagggcGAGCAGAGAGCAGAGCATCC-3' 
with T3 and T7 overhangs for sequencing (in lower case). PCR 
reactions (0.5 x HF buffer, 0.16 oM dNTPs, 0.5 uM primers, 
0.8 U Phusion Hi-Fi polymerase (Thermo Scientific), 100 ng 
DNA) were incubated in the PTC DNA Engine thermocycler 
(Bio-Rad, Hercules, CA, USA) under the following conditions: 
95°C for 3 min; 30 cycles of 95°C for 20 s, 60°C for 30 s, 72°C 
for 1 min; 72°C for 5 min. PCR fragments were purified using 
the QIAquick PCR purification kit following the manufacturer's 
recommendations and Sanger-sequenced. The sequences were 
assembled using SeqMan software (DNASTAR Lasergene® 8 
Core Suite, Madison, WI, USA). 

Multiple alignment and phylogeny reconstruction 

A barley cDNA, HvLUXl homologous to HvLUXl was extracted 
from NCBI GenBank using Blastn search in the nr/nt database 
and, because of the misannotation of the LUX2 cDNA in Gen- 
Bank, it was manually translated. HvLUXl and HvLUX2 pro- 
teins were further used as NCBI Blastp queries to extract 
homologs from other plant species (E-value cutoff 10 -30 ; last 
search on 22 November 2012). In addition, LUX-like protein 
from S. bicolor (Sb03 g039610) was retrieved from the Gramene 
database (http://www.gramene.org/Sorghum_bicolor). One rep- 
resentative per taxa was kept in case several identical sequences 
were retrieved. The nomenclature and accession numbers of ana- 
lyzed LUX-like sequences (species - genotype (database accession 
number)) are as follows: A. thaliana: Atl — Columbia (LUX1, 
NP_001 190022; NOXla, NP_200765), At2 - Wassilewskija/ 
Landsberg erecta (NOXlb, AAM65635); A. lyrata: Al - MN47 
(LUX1, XP_002875796; NOX1, XP_002864640); B. 
distachyon: Bd - Bd21 (LUX1, XP_003565l6l; LUX2, 
XP_003 567265); Glycine max. Gm - Williams 82 (LUX1, 
XP_003539607; LUX2, XP_003 537977); Lycoris longituba: LI - 
six mixed varieties (LUX1, ADG58105; LUX2, ADG58070; 
LUX3, ADG57915; LUX4, ADG57989); Medicago truncatula: 
Mt - A17 (LUX1, XP_003606714); Nicotiana benthamiana: Nb 
— unknown (LUX1, BAE16280); H. vulgare ssp. vulgare. Haruna 

New Phytologist (2013) 199: 1045-1059 
www.newphytologist.com 



048 Research 



New 

Phytologist 



Nijo (LUX1, BAJ88719; LUX2, AK356714 manually trans- 
lated); Hyacinthus orientalis: Ho — Delft Blue (LUX1, 
AAS21003); Oryza sativa ssp. indica: Osi — 93-11 (LUX1, 
EAY77400; LUX2, EEC71789); O. sativa ssp. japonica: Osjl - 
Zhonghual 1 (LUXla, AAS90600), Osj2 - Nipponbare (LUXlb, 
NP_001045537; LUX2, NP_001044783); Physcomitrella patens 
ssp. patens. Pp - Gransden2004 (LUX1, XP_00 1757337; 
LUX2, XP_001768612; LUX3, XP_001757676; LUX4, 
XP_00 1757677); Populus tricbocarpa: Pt - Nisqually-1 (LUX1, 
XP_002299857; LUX2, XP_0023 14163); Ricinus communis-. Rc 
- Hale (LUX1, XP_002520534); Selaginella moellendorffii: Sm - 
Plants Delights Nursery (LUX1; XP_002962201); Solanum 
lycopersicum: SI — unknown (LUX1, BAE16281); S. tuberosum: 
St - unknown (LUX1, BAE16282); S. bicolor. Sb - BTx623 
(LUX1, XP_002459177; LUX2, Sb03 g0396l0); Vitis vinifera: 
Vv - PN40024 (LUX1, XP_002283159); Zea mays: Zm - mixed 
genotype (LUXla, NP_001143908; LUXlb, NP_001 147359; 
LUX2, DAA56868). Protein sequences were aligned using 
MAFFT v.6.851b with the E-INS-i algorithm (Katoh & Toh, 
2008). Visually unaligned regions were cropped from the 
alignment using BioEdit 7.0.9.0 software (Hall, 1999). The 
alignment was converted from fasta to Phylip4 format using 
ReadSeq 2.1.30 tool implemented at EBI (http://www.ebi.ac.uk/ 
cgi-bin/readseq.cgi). A best-scoring maximum likelihood tree and 
pairwise distances were calculated using the fast bootstrapping 
algorithm of the RAxML 7.3.5 with the GAMMAJTT model 
and 'autoMRE' parameter determining an optimal number of 
bootstrap replicates (Stamatakis etal, 2008). The best-fitting 
amino-acid substitution model (JTT) was estimated using Pro- 
teinModelSelection Perl script (http://www.exelixis-lab.org/soft- 
ware/ProteinModelSelection.pl). The tree was edited in 
Dendroscope 3 software (Huson & Scornavacca, 2012). Micro- 
synteny analysis at the LUX loci in B. distachyon, O. sativa ssp. 
japonica, and S. bicolor was performed using CoGe GEvo script 
(Lyons & Freeling, 2008) with the default BlastZ parameters 
except for the score threshold 6000 (Methods S4). The sequence 
logo of the SHAQKYF motif was generated using WebLogo 
(Crooks etal., 2004). 

Analysis of HvLUXI natural diversity 

A CDS of HvLUXI was amplified and sequenced as indicated in 
the 'Segregation analysis' section (GenBank accession nos. 
KC668259-KC668274). The sequences were assembled using 
SeqMan software (DNASTAR Lasergene* 8 Core Suite) and 
polymorphic haplotypes (sequence variants) were discerned man- 
ually based on the revealed polymorphisms. The Median Joining 
haplotype network was constructed using SplitsTree4 with 
default parameters (Huson & Bryant, 2006). Indels were 
encoded as a nucleotide substitution to include them in the anal- 
ysis as an informative character. 

Statistical analysis 

Significant differences in flowering time among Bowman, 
V>owmzn{eamlO), and Bowman(^wi0+ Ppd-Hl) grown under 
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SDs or LDs were calculated using paired £-test with a 95% 
confidence level (P<0.05). Significant differences in meristem 
development and gene expression between Bowman and 
Bowman(^wiO) were calculated using a general linear model 
in the SAS software 9.1.3 (SAS Institute, 2009) with the 
factors genotype, time point, biological replicate, and first- 
order interaction effects. Significant differences (P<0.05) 
between least-squares means (LSmeans) of the genotype x time 
interactions were calculated using a Tukey-Kramer adjustment 
for multiple comparisons. 

Results 

The eamW mutant is early flowering under SD and LD con- 
ditions 

In this study, we analyzed the effect of earn 10 on photope- 
riod-dependent flowering in barley and its interaction with 
natural variation in the major barley photoperiod response 
gene Ppd-Hl. Flowering time was scored in the spring cv 
Bowman(ppd-Hl) and in three derived introgression lines 
Bowman(eamlO), Bowman(Ppd-Hl), and Bowman(/^- 
H1+ eamlO) under LDs and SDs (Fig. la). Under LDs, Bow- 
mzn{Ppd-Hl+ eamlO) flowered first at 27 d after sowing 
(DAS), followed by ¥>awman.{Ppd-Hl) at 31 DAS, Bowman 
(eamlO) at 42 DAS, and Bowman at 46 DAS. Under SDs, 
Bowman and Bowman(Ppd-Hl) did not flower until 90 DAS 
when the experiment was terminated. By contrast, Bowman 
(eamlO) and Bowman(Ppd-Hl+eamlO) flowered under SDs at 
76 and 35 DAS, respectively. Thus, earn 10 accelerated flower- 
ing under both LDs and SDs. The combination of Ppd-Hl 
and earn 10 resulted in the earliest flowering phenotype irre- 
spective of the photoperiod. Interestingly, variation at Ppd-Hl 
caused a significant difference in flowering time in the back- 
ground of eamlO under SDs. Taken together, these findings 
suggest that Ppd-Hl interacts epistatically with earn 10 to accel- 
erate flowering under both LD and SD conditions. 

To determine the effect of eamlO on the development of the 
shoot apical meristem (SAM), we dissected meristems of 
Bowman and Bowman(&zm70) plants starting 1 wk after germi- 
nation until heading and scored morphological changes of the 
SAM based on the Waddington developmental scale (Fig. 2a; 
Waddington etal, 1983). Under LDs and SDs, the SAM of 
Bowman(eamlO) developed significantly faster than that of Bow- 
man (Table SI). Significant genetic differences at single time 
points of dissection were only detected after Waddington stages 8 
and 4 in LD and SD conditions, respectively (Fig. 2a). Under 
SDs, the most prominent differences between genotypes were 
observed in the distance between the base and the meristem (stem 
elongation), which was significantly larger in BowmsLn(eamlO) 
than in Bowman at most time points (Fig. 2a). Therefore, eamlO 
primarily affected stem elongation and had marginal effects on 
SAM development. 

In barley, expression of HvFTl correlates with flowering 
time (Turner etal, 2005). To test whether differences in the 
expression of HvFTl could explain the observed flowering-time 
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Fig. 1 Flowering time and HvFT1 expression levels of barley (Hordeum vulgare) wildtype and introgression lines. Bowman and introgression lines carrying 
PHOTOPERIOD 1 gene (Ppd-H1), early maturity 10 (eamW) or Ppd-HI + eamW grown under long days (open bars) or short days (closed bars), (a) 
Flowering time is shown as mean days to emergence of the main spike awns. Error bars represent SD for 15-18 plants of each genotype. Genotypes with 
the same letter are not significantly different by paired f-test. (b) HvFT1 expression relative to HvActin is shown as average of two biological and two 
technical replicates ± SD. Leaf material was sampled after 2 wk (long days) and 4 wk (short days) 2 h before light-off. ND, not detectable. 



phenotypes, we measured expression of HvFTl in the different 
Bowman introgression lines under LDs and SDs (Fig. lb). 
After 2 wk under LDs, HvFTl was expressed only in lines with 
the dominant Ppd-HI allele, Bowman(Ppd-Hl) and Bowman 
(Ppd-HI +eam 10). After 4wk under SDs, HvFTl expression 
was detected in Bowman(Ppd-Hl +eaml0), but not in the 
other genotypes. To identify expression differences between 
Bowman and Bowman (earn 10) during plant development, we 
analyzed HvFTl expression in leaf samples from plants dis- 
sected for meristem scoring under LDs and SDs (Fig. 2b). 
Under LDs, HvFTl expression was higher in Bowman(<M»27(9) 
than in Bowman at the double-ridge stage until early stem 
elongation (Waddington stages 2—6; 9 and 24 DAS). Under 
SDs, HvFTl expression was detected in Bowman(^w70) after 
the beginning of stem elongation (beyond Waddington stage 3; 
79 DAS), but not in Bowman. Thus, the eamW mutation 
caused an up-regulation of HvFTl and accelerated flowering 
time under both LDs and SDs. The dominant Ppd-HI allele 
in the background of eamlO caused a stronger up-regulation of 
HvFTl than the recessive ppd-hl allele under both LD and 
SD conditions (Fig. lb). 

eamlO affects the expression of circadian clock homologs 

Previous studies have shown that mutations in circadian clock 
genes caused earliness and day-neutrality in Arabidopsis and bar- 
ley plants (Davis, 2002; Hazen etal, 2005; Onai & Ishiura, 
2005; Faure etal., 2012; Zakhrabekova etal., 2012). To investi- 
gate whether the eamlO mutation led to a disruption of clock 
rhythmicity, we studied diurnal (under SDs) and circadian (LL) 
expression of barley clock homologs, HvCCAl, Ppd-HI 
{HvPRR37X HvPRR73, HvPRR59, HvPRR95, HvPRRl, HvGI, 
and HvLUXl in Bowman and Bowman(m»27(9) (Fig. 3). 

Under SDs, Bowrmn(eaml 0) showed reduced expression of 
HvCCAl and HvPRR73 compared with Bowman. Conversely, 
expression levels of Ppd-HI and HvPRRl were higher in Bow- 
miLn(eamlO) than in Bowman during the day, while the night 
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peak of expression of HvLUXl and HvPRR73 was higher in 
Bowman(^w70) than in Bowman. Expression differences 
observed under SDs were even more pronounced under LL. 
The strongest differences between Bowman and the eamlO 
mutant were detected in the expression of HvCCAl and Ppd- 
HI under LL. Under these conditions, in Bowman(<M»27(9), 
circadian amplitude of HvCCAl expression was significantly 
reduced during subjective days, while Ppd-HI was significantly 
up-regulated in Bowman(crfw7(9) compared with Bowman at 
all time points. In contrast to Ppd-HI, expression of its homo- 
logs HvPRR73, HvPRR59, and HvPRR95 was not strongly 
affected by the eamlO mutation. The evening-expressed genes 
HvPRRl and HvGI were up-regulated in Bowman(m»27(9) 
compared with Bowman. Taken together, these results indicate 
that eamlO alters the expression of barley homologs of Arabid- 
opsis clock genes, in particular the expression of Ppd-HI and 
HvCCAl. 

In Arabidopsis, the expression of genes implicated in the regu- 
lation of photosynthesis and photoperiod-dependent flowering is 
under circadian control. We therefore tested whether eamlO 
changed the diurnal and circadian expression of the following 
genes: the photoperiod response genes HvCOl and HvFTl; the 
chlorophyll A/B binding protein gene HvCABIIl from the photo- 
synthetic pathway; and HvCCR2, encoding the barley ortholog 
of the GYLCINE-RICH RNA-BINDING PROTEIN 7 (GRP7I 
CCR2), characterized as a slave (nonself-sustaining) oscillator 
(Scheming & Staiger, 2005; Fig. SI). Under SDs, HvCOl did 
not show differences in expression, and HvFTl expression was 
not detected in any of the two genotypes at 21 DAS. Under LL, 
HvCOl showed a small but significant up-regulation in Bowman 
(eamlO) at some time points during the subjective night. Bow- 
man(eamlO) exhibited a strong up-regulation of HvFTl during 
the subjective day compared with Bowman. In Bowman(eamlO), 
HvCABIIIwus up-regulated at nights and subjective nights, while 
the peak expression of HvCCRl was reduced under SDs. The 
earn 10 mutation thus affected the expression of clock output 
genes. 
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Fig. 2 Meristem development and expression levels of HvFTI under long days and short days in barley (Hordeum vulgare) Bowman and Bowman(eam70) 
plants, (a) Apex development assessed by the Waddington scale under long- and short-day conditions in Bowman (black lines) and Bowman(eam70) (red 
lines). Timings of double ridge formation (WS2) and the initiation of internode elongation (WS3) are indicated by horizontal full and dashed lines, 
respectively. Black and red arrows indicate the day when Bowman and Bowman(eam70) headed under long-day conditions. The small chart above shows 
the distance in cm between the base of the stem and the base of the meristem measured under short-day conditions on the indicated days. Values are 
means of three plants ± SD. Significant differences in meristem development and stem elongation are indicated: *, P<0.05. (b) HvFTI expression relative 
to HvActin in Bowman (black lines) and Bowman(eam7Q) plants (red lines) under long- and short-day conditions. The second last leaves of three plants 
were harvested and pooled. Sampling was done 2 h before light-off. The average Waddington stage of meristems at the respective sampling times is 
reported in the table. Values represent the average of two biological and two technical replicates ± SD. 



Identification of barley Hvluxl as a candidate gene under- 
lying eamW 

To identify candidate genes underlying the earn 10 QTL, we 
investigated polymorphism of 133 flowering-related barley genes 
in Bowman and Bowman(mw70) using high-throughput 
sequencing of libraries enriched for flowering-related genes. In 
total, we found five SNPs between Bowman and Bowman 
(eamlO) (Fig. S2). These SNPs resided within three flowering- 



related HarvEST unigenes, 16001, 19636, and 22370, homolo- 
gous to the Arabidopsis genes CONSTITUTIVE PHOTO- 
MORPHOGENIC 1 (COP1), LUX, and ABA-RESPONSIVE 
ELEMENT BINDING PROTEIN 2 (ABE4), respectively. 

As reported by Druka etal. (2011), Bowman(ft?w70) carries 
introgressions on four chromosomes, 3H (single polymorphic 
marker), 4H (c. 1 1 cM), 6H (c. 15 cM), and 7H (single polymor- 
phic marker). Using the GenomeZipper (Mayer etal., 2009), we 
mapped all three barley genes polymorphic between Bowman 
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Fig. 3 Expression patterns of main circadian clock genes in barley (Hordeum vulgare). Expression of HvCCAl , HvLUXl, Ppd-Hl (HvPRR37), HvPRR73, 
HvPRR59, HvPRR95, HvPRRI , and HvGI in Bowman (black lines) and Bowman(eam70) plants (red lines) under short-day and continuous-light conditions. 
White, black and grey bars indicate days, nights and subjective nights, respectively. Values represent averages of two biological and two technical 
replicates of expression values relative to HvActin ± SD. Significant differences in gene expression are indicated: *, P < 0.05. 



and Howman(eamlO) within the introgressed regions. The barley 
LUX-like gene designated as HvLUXl was located on the linkage 
group 3H (Hv3) at the position of POPA marker 2_1500 
(171.64 cM), whereas the two other barley genes were mapped 
on the linkage group 6H both in the region from 85.16 to 
91.79 cM. Among the three genes, only HvLUXl colocated with 
the earn 10 QTL, mapped on the distal end of the long arm of 
chromosome 3H below the marker ABC166 (155.85 cM; Borner 
etal, 2002) (Fig. S2). 

The Harvest unigene 19636 comprised only a partial coding 
sequence of HvLUXl. Therefore, we extracted its homolog 
Hv.20312, a full-length unigene, from NCBI GenBank 
(AK357505). Similarly to the Arabidopsis LUX gene (Hazen 
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etal, 2005; Onai & Ishiura, 2005) characterized as a transcrip- 
tion factor (TF), the intronless HvLUXl gene encodes an MYB- 
domain-containing SHAQKYF-type GARP family protein of 
274 amino-acid residues (Fig. 4a). An evening element (EE, 
AAATATCT) characteristic of clock genes expressed in the even- 
ing (e.g. Arabidopsis TOC1 and LUX, Alabadi etal., 2001; 
Hazen etal, 2005) was found 208 bp upstream of the start 
codon. In addition to the EE, the promoter region of HvLUXl 
comprised two conserved LUX binding site motifs (GATACG 
and GATTCG) required for L UX autoregulation in Arabidopsis 
(Heifer etal., 2011). 

Database search identified 36 LUX-like proteins from 19 plant 
species. Multiple alignment of these proteins revealed two highly 
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conserved regions corresponding to the MYB and as yet unchar- 
acterized domains (Fig. S3). The family-defining SHAQKYF 
motif (SHLQKY(R/Q) in the case of Z£ZY-like genes) was con- 
served across all LUX-like sequences, thus indicating that the 
integrity of this motif may be critical for the function of LUX- 
like TFs. The mutation identified in the LUX-like gene from 
Howman(eamlO) changed the first amino-acid residue in this 
conserved motif (Fig. 4b; GenBank accession no. KC668258). 
Therefore, we suggested that this could lead to a disruption of 
LUX function and tentatively assigned Hvluxl as a candidate 
gene underlying earn 10. 

To corroborate the link between the Hvluxl mutation and the 
early flowering, we conducted a segregation analysis using 1002 
BCiF 2 lines derived from backcrossing Bowman(<M»z70) to Bow- 
man. Based on sequencing analysis of early-flowering BCiF 2: 3 
lines, we confirmed that the Hvluxl mutation cosegregated with 
the earliness phenotype. 

Since the region of the eaml 0 QTL, spanning c. 17 cM, might 
contain alternative candidate genes, explaining the observed phe- 
notypes and physically linked to Hvluxl, we searched this region 
for homologs of clock and flowering genes. A barley homolog of 
Arabidopsis response regulator {ARB) clock gene (Salome et ai, 
2006) was found c. 16cM above HvLUXl and excluded as an 
alternative candidate based on the sequence analysis (Notes SI; 
Tables S4, S5). 

Natural diversity of HvLUXl 

Natural variation at flowering time loci in barley has been pro- 
posed to be adaptive (Jones et ai, 2008; Comadran etal, 2012). 
To investigate natural variation in HvLUXl, we sequenced the 
full-length gene from 88 diverse wild and cultivated barley 
genotypes from various geographic origins. A haplotype analysis 



identified 16 HvLUXl haplo types (Fig. 5a) distinguished by 
seven SNPs and six indels. The nonsynonymous SNP(545) and 
six in-frame indels resulted in minor changes in the amino-acid 
composition of HvLUXl, although none of them led to a protein 
truncation. The majority of the investigated genotypes, 52 sam- 
ples, carried a 'haplotype 1' of HvLUXl. Diversity of HvLUXl 
was much higher in wild than in cultivated genotypes. We identi- 
fied 16 haplotypes in 52 wild barley accessions, while 35 out of 
36 cultivated accessions comprised a 'haplotype 1'. In our study, 
mutations in HvLUXl did not associate with particular subspe- 
cies or place of origin (Fig. 5b; Table S2). The Hvluxl allele was 
not found in the natural population, apparently because of its 
artificial origin, but the pattern of its variation clearly suggests 
that it originated from haplotype 1 . 

Phylogenetic framework of LUX-\\ke genes in plants 

The maximum-likelihood analysis of the LUX homologs 
extracted by the Blast search defined a LUX-family clade com- 
prising Arabidopsis LUX along with the 38 LUX-like proteins 
from 7 monocots, 10 dicots, the spikemoss S. moellendorffii, and 
the moss P. patens. To infer phylogenetic relationships within the 
LUX family, we recalculated the phylogeny, including only mem- 
bers of the LUX-family clade (Fig. 6). The LUX-like proteins 
from the angiosperms fell into two monophyletic clades of dicots 
and monocots. Within these clades, distinct subclades and 
branches at the level of orders and families could be recognized. 
The topology of the Brassicales and Poales clades was indicative 
of independent duplications of the ZC/X-like loci in these taxa. In 
the Arabidopsis genus, such duplication apparently preceded spe- 
ciation {LUX and NOX/BROTHER OF LUX genes; Hazen etal, 
2005; Dai etal., 2011), whereas in monocots we tentatively 
ascribed it to a deeper phylogenetic level, the common ancestor 
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Fig. 4 The structure of the barley (Hordeum vulgare) HvLUXl gene and the effect of a point mutation, (a) The structure of the HvLUXl gene and the 
position of a nonsynonymous single nucleotide polymorphism (SNP) in the MYB domain. The exon region is shown in blue with the region encoding the 
GARP family MYB domain (approximate location) highlighted in red, and the UTR or promoter regions shown in grey. PCR primers for amplification of 
HvLUXl are shown as horizontal arrows. Distances are shown in pairs of nucleotides relative to the full-length HvLUXl gene (NCBI accession no. 
AK357505). A vertical arrow indicates the position of a nonsynonymous SNP with the wild allele (A) in HvLUXl and the mutated allele (T) in Hvluxl. An 
evening element (*, AAATATCT) and two LUX binding site motifs (**, GATACG; ***, GATTCG) are shown in the promoter region, (b) The amino-acid 
mutation in the highly conserved SHLQKY(R/Q) motif in Hvluxl . An alignment of 37 SHLQKY(R/Q) motifs from 19 species (Fig. S3) is showed as a 
sequence logo. Green letters, polar amino-acid residues (aa); blue, basic aa; red, acidic aa; black, hydrophobic aa. The height of the letters illustrates the 
frequency of a corresponding aa in the alignment. 
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Fig. 5 HvLUXl haplotypes in wild and cultivated barley (Hordeum vulgare). (a) Sixteen polymorphic haplotypes of HvLUXl. Insertions: INS1 , 
GGAGGAGGA; INS2, ACAGCAACA; INS3, ACA; INS4, ACAGCAACAACA; INS5, GCG; deletion: DEL1 , GTC. A nonsynonymous single nucleotide 
polymorphism (SNP) is highlighted in red. (b) Median-joining network of 16 HvLUXl haplotypes. Numbers at the nodes indicate the number of genotypes 
carrying the corresponding haplotype (out of 88 accessions). The haplotype frequency is also reflected in the relative size of a node. The color of a node 
corresponds to the different species: blue, Hordeum vulgare ssp. vulgare; red, H. vulgare ssp. spontaneum; green, H. agriocrithon. The mutated X-ray- 
derived Hvluxl allele from Bowman is labeled with an asterisk. 



of the Poaceae species (duplicates assigned as LUX1 and LUX2). 
It is remarkable that the LUX outparalogs {sensu Sonnhammer & 
Koonin, 2002) in Poaceae significantly diverged after duplication 
as manifested by the distance between the LUX1 and LUX2 
clades (ML distance 1.22). By contrast, the divergence between 
the Arabidopsis LUX-like outparalogs remained relatively low 
(ML distance 0.42). 

To investigate in detail the ancestral LUX duplication in the 
Poaceae, we analyzed microsynteny of the regions comprising 
LUX-Yike loci in Brachypodium, barley, rice, and sorghum 
(Fig. 7). This analysis revealed a 'fuzzy' homology of a c. 40 kbp 
fragment of the barley physical map (IBSC, 2012), carrying the 
HvLUXl locus, with the corresponding genomic regions of the 
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other Poaceae species. Apparently, the order of the contigs in this 
region of the barley physical map awaits further refinement. Nev- 
ertheless, we could locate HvLUXl on the barley chromosome 
3H (Hv3). HvLUXl was not found on either the physical or the 
genetic maps of barley. 

The LUX1 genes resided in syntenic regions of the chromo- 
somes Bd2, Hv3, Osl, and Sb3 of Brachypodium, barley, rice, 
and sorghum, respectively, which have been shown to be ortholo- 
gous (Mayer etal, 2011). It is noteworthy that a c. 10 kbp frag- 
ment comprising BdLUXl was inverted, while the immediately 
adjacent orthologous gene pairs were in a collinear order. The 
LUX2 gene was found on the same chromosomes in the inverted 
fragment located c. 5.2—6.8 mbp upstream of the LUX1 locus in 

New Phytologist (2013) 199: 1045-1059 
www.newphytologist.com 



054 Research 



New 

Phytologist 



4g j-Osj2LUXlb 
JOsiLUXl 



Osj2LUX2 




Asparagales 

Malpighialis 

Fabales 

Vitales 

Solanales 

Brassicales 



-PpLUXl 



Fig. 6 Maximum-likelihood phytogeny of LUX-like proteins in 1 9 plant species. Bootstrap values (%) are shown at the nodes. Clock pictograms mark 
Arabidopsis paralogs, LUX and NOX implicated in regulation of the circadian clock (Dai etal., 2011; Heifer eta/., 2011). The LUX-like protein encoded by 
the gene described in this study is marked by a red star pictogram. Scale bar, 0.1 amino-acid substitutions per site. Colour of the branches indicates 
angiosperm orders. Protein nomenclature is as described in the Materials and Methods section. 



all the three grass species. Several genes in the vicinity of LUXloci 
were unique to either rice or Brachypodium. This implies that 
gene gain and loss events occurred at these loci after speciation. 

Discussion 

Bowman(eam70) is a circadian clock mutant 

In this study, we describe the barley mutant locus earn 10, which 
accelerates flowering under LDs and, in addition, allows the plant 
to flower and mature under SDs. eamlO has been isolated and 
described by Favret & Ryan (1966), when early cultivars had 
become an important goal in barley breeding programs and a 
large number of additional earn mutants were isolated in the 
Swedish breeding programs (Lundqvist, 2009). Recently, muta- 
tions underlying eam.8 have been located in HvELF3 (Faure etal., 
2012; Zakhrabekova etal., 2012). These mutations resulted in 
the disrupted expression rhythmicity of barley homologs of Ara- 
bidopsis clock genes such as HvCCAl, HvPRRl and HvGI, the 
up-regulation of genes involved in the photoperiodic pathway 
(e.g. Ppd-Hl) and the acceleration of flowering independent of 
photoperiod. Similar effects of mutations in ELF3-\ike genes on 
the rhythmicity of clock genes and photoperiodic reactions have 
been reported in Arabidopsis, rice, pea (Pisum sativum), and 
lentils {Lens culinaris) (Zagotta etal., 1996; Saito etal., 2012; 
Undurraga etal., 2012; Weller etal, 2012; Zhao etal, 2012). In 
Bowman(wwIO), diurnal and circadian expression patterns of 
clock components such as HvCCAl, HvPRR95, and HvGI 
dampened under LL conditions, thus indicating that the 
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circadian rhythms in Bowman(^w70) were compromised. In 
both <Mw8and eamlO mutants, disruption of circadian rhythmic- 
ity, as attested by the change in expression patterns of clock 
genes, accompanied acceleration of flowering irrespective of day 
length. As shown for Bowman(fiiwSJ (Faure etal., 2012), early 
flowering of 'Qowman{eaml0) was associated with an upregula- 
tion of HvFTl under SDs (Fig. 2b). Thus, mutations in eam8 
and eamlO caused an induction of the long-day photoperiod 
pathway under noninductive SD conditions. Therefore, eam8 
and earn 10 might be part of the mechanism modulating light sig- 
nal transduction from receptors to downstream components of 
the photoperiod pathway as has been shown for the EC genes in 
Arabidopsis (McWatters etal, 2000; Carre, 2002). 

Rate of stem elongation was significantly accelerated in earn 10 
mutants, whereas the timing of meristem development closely 
followed that of Bowman, especially at the earlier stages. In Ara- 
bidopsis, abnormally elongated hypocotyls are characteristic of a 
dysfunctional circadian clock, and molecular pathways underly- 
ing this phenomenon have been outlined (Niwa etal, 2009; 
Nusinow etal, 2011). It is tempting to speculate that while 
organs abnormally elongated at the early developmental stages of 
circadian-defective barley and Arabidopsis plants are clearly dif- 
ferent, underlying molecular mechanisms might be similar. 

eamlO regulates expression of Ppd-H1 

Measurements of circadian-regulated gene expression suggested 
that the circadian defects in Bov/man(eaml 0) were weak under 
SD conditions and significantly enhanced under LL in the 
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Fig. 7 Microcollinearity in the vicinity of duplicated LUX-Wke loci in Poaceae species. Comparison of the Sorghum bicolor chromosome 3 (Sb3), the Oryza 
sativa ssp. japonica chromosome 1 (Os1), and the Brachypodium distachyon chromosome 2 (Bd2) was based on the microsynteny between corresponding 
regions on the physical maps of these species (Sb, JGI map v1 .4; Os, MSU Rice Genome Annotation v7; Bd, JGI v1) as revealed by the CoGe GEvo script 
(http://genomevolution.org/CoGe/GEvo.pl). Location of the Hordeum vulgare ssp. vulgare LUX1 locus (HvLUXl, as POPA marker 2_1 500) is shown on 
the consensus genetic map (as found in the barley GenomeZipper); location of HvLUX2 could not be established. Orthologs are connected by horizontal 
solid lines, outparalogs by dashed lines. LUX-Wke genes (SbLUXI, Sb03 g047330; Osj2LUX1b, LOC_Os01 g7420; BdLUXl , Bradi2 g62070; SbLUX2, 
Sb03g039610; Osj2LUX2, LOC_Os01 g62660; BdLUX2, Bradi2 g54790) are highlighted in red. Scale bar, 10kbp. 



absence of external cues. Light and temperature signals could thus 
drive rhythms in ^owmnn(eamlO). However, under LD and SD 
conditions when external cues were present, eamlO had 
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significant effects on flowering time, apparently through a direct 
or indirect interaction with Ppd-Hl. eamlO mutant lines with 
the wildtype Ppd-Hl allele flowered earlier than those with the 
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recessive ppd-Hl allele under LDs and SDs. This difference was 
most pronounced under SDs. Thus, variation at Ppd-Hl in the 
background of eamlO, affected flowering under SDs, while previ- 
ously it has been postulated that Ppd-Hl acts exclusively under 
LDs (Turner etai, 2005). The interaction between eamlO and 
Ppd-Hl was further supported by a significant up-regulation of 
Ppd-Hl expression in Bowman(«zm70) compared with Bowman 
under both SDs and LLs. These results suggested that EamlO acts 
as a repressor of Ppd-Hl. Early flowering of lines carrying Ppd- 
Hl and eamlO correlated with an up-regulation of HvFTl under 
SDs (Fig. lb), presumably caused by the up-regulation or differ- 
ences in the diurnal expression pattern of Ppd-Hl in the eamlO 
background. 

HvLUXI is a candidate gene underlying eamlO 

The identification of a mutation in an extremely conserved 
region of HvLUXI strongly suggested Hvluxl as a gene underly- 
ing the eamlO locus. The segregation analysis, the flowering and 
expression phenotypes, and the screening of genes in the vicinity 
of earn 10 solidified the link between Hvluxl and earn 10, despite 
the possible shortcomings of using introgression lines carrying 
multiple genes transferred from donor germplasm. The mecha- 
nistic role of Hvluxl in the eamlO flowering and expression phe- 
notypes is yet to be discovered. 

It has been shown in Arabidopsis that LUX and ELF3 together 
with ELF4 form the so-called 'evening' complex and have the 
same transcriptional targets (Kolmos etal, 2011; Nusinow etal, 
201 1; Herrero etal, 2012). Both LUX and ELF3 physically asso- 
ciate with the promoter of PRR9 to repress its transcription 
(Dixon etal, 2011; Herrero etal, 2012). Up-regulation of Ppd- 
Hl in eamlO and eam8 mutants (Faure etal, 2012) suggested 
that HvLUXI and HvELF3 act in the same pathway, and that 
their transcriptional targets are partly conserved between Arabid- 
opsis and barley. Interestingly, the repressive function of eamlO 
was only seen for HvPRR37, while expression of its paralogs 
HvPRR73 and HvPRR59/95 was not strongly affected. By con- 
trast, EC genes in Arabidopsis control expression of PRR7 and 
PRR9 (Kolmos etal, 2009; Herrero etal, 2012). Differences in 
the transcriptional targets of barley Eam8/Eaml0 and Arabidopsis 
ELF3 suggest a different clock construction in barley and Arabid- 
opsis. Bov/man(eamlO) was characterized by a down- regulation 
of HvCCAl and an up-regulation of HvPRRl, HvGI, and 
HvLUX, This is consistent with the role of Arabidopsis LUX, 
which activates CCA1 and represses TOC1 (PRR1), GI, and 
LUX itself (Hazen etal, 2005; Onai & Ishiura, 2005). Up-regu- 
lation of HvLUX m earn 10 and the presence of two conserved 
LUX binding site motifs (GATACG and GATTCG) in the pro- 
moter region of HvLUX suggested that HvLUXI controls its 
own expression, as seen in Arabidopsis (Heifer etal, 201 1). 

LUX duplicated independently in Poaceae and Arabidopsis 

Interspecific comparison of clock genes within phylogenetic 
frameworks may facilitate integration of knowledge and also pin- 
point dissimilar patterns of gene evolution, when direct 
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comparisons between homologous genes from distant plant lin- 
eages should be made with greater caution. Phylogeny-supported 
comparison of flowering genes in grasses and Arabidopsis has 
identified several cases when direct orthology between members 
of gene families could not be established (e.g. CDF1, TOE1, 
TEM, HAP genes; Higgins etal, 2010). In the case of clock 
genes, such information is extremely scarce. It has been suggested 
that, in monocots and eudicots, genome duplication events led to 
the independent duplication and diversification of the PRR gene 
family (Murakami etal, 2007; Takata etal, 2010; Campoli 
etal, 2012b). 

In our study, LUX-Xihz genes could be traced down to mosses 
and apparently represent an ancient branch of the GARP tran- 
scription factor family. We found duplication of LUX genes in 
the ancestor of Poaceae (LUX1 and LUX2) independent of the 
duplication in Arabidospis (LUXand NOX). While ZMfparalogs 
in Arabidopsis retained partially redundant functions (Dai etal, 
2011) and remained phylogenetically close, LUX paralogs in 
grasses significantly diverged before grass species radiation. This 
suggests that evolutionary fate of individual members of clock 
gene families might be different, especially in distant plant 
lineages. 

Microsynteny around LUX loci in Brachypodium, rice, and 
sorghum revealed ancestral duplication, which is consistent with 
the phylogeny. However, we could not resolve whether this 
duplication originated from the ancient whole-genome duplica- 
tion event or local translocation. Several gene loss/gain and micro- 
inversion events occurred independently in the genomes of the 
three grass species, while both LUX \oa remained intact. These 
data implicate selection for retention of LUX loci in grass 
genomes, corroborating previous findings, which suggested pref- 
erential retention of duplicated clock genes in Brassica rapa (Lou 
etal, 2012). 

HvLUXI is under purifying selection 

Structured allelic variation at several flowering- and clock-related 
loci has been proposed as a signature of plant adaptation to either 
natural or agricultural environments (Izawa, 2007; Clotault etal, 
2012; Weller etal, 2012; Anwer & Davis, 2013). Allelic varia- 
tion at three flowering loci, Ppd-Hl, HvCEN, and HvELF3, asso- 
ciated with the geographical distribution of wild and cultivated 
barley habitats has been suggested to be adaptive (Jones etal, 
2008; Comadran etal, 2012; Faure etal, 2012; Zakhrabekova 
etal, 2012). 

We found that HvLUXI diversity is high (16 haplo types in 88 
accessions). Surprisingly, only one haplotype was prevalent in 
cultivars and landraces from diverse geographic origins; the diver- 
sity of HvLUXI was predominantly confined to wild barley (16 
haplotypes in 49 accessions). Remarkably, all in-frame indels 
identified in several HvLUXI haplotypes resulted from variation 
in a number of microsatellite repeats. Functional and adaptive 
variation in the number of tandem amino-acid repeats have been 
demonstrated in flowering and clock genes, such as ELF3 in Ara- 
bidopsis and FCA in endemic Hawaiian mint {Stenogyne rugosd) 
(Lindqvist etal, 2007; Undurraga etal, 2012). Whether 
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variation in HvLUXl is functional warrants further investigation. 
The absence of frameshifts and predominance of a single haplo- 
type indicate that HvLUXl was under purifying selection. While 
barley plants with nonfunctional HvelfB might have selective 
advantage in certain environments, integrity of HvLUXl seems 
to be under functional constraints. 

Conclusion 

Isolation and characterization of Hvluxl advance our understand- 
ing of the photoperiod response pathway in temperate cereals. 
The mutation in HvLUXl is linked to eamlO, which caused early 
flowering independent of the photoperiod. Hvluxl was associated 
with similar, but weaker, molecular and flowering phenotypes 
than HvelfB. It thus represents an interesting resource for modu- 
lating photoperiod response and ultimately adaptation in barley. 
Revealed patterns of diversification of duplicated LUX loci and 
their natural variation offer some insight into a number of ques- 
tions regarding the differential fate of duplicated genes in differ- 
ent plant lineages and the impact of the circadian clock 
components on adaptation in natural populations. 

Acknowledgements 

M.v.K. and S.J.D. are supported by the Max Planck Society and 
by grants from the DFG SPP1530 ('Flowering time control: 
from natural variation to crop improvement'). A.P. was sup- 
ported by the International Max Planck Research School fellow- 
ship. We are grateful for seeds of barley lines from Benjamin 
Kilian (IPK-Gatersleben) and Eyal Fridman (Hebrew University 
of Jerusalem). We cordially thank Kerstin Luxa, Teresa Bisdorf 
and Elisabeth Luley for excellent technical assistance and Ganga 
Jeena for help with the data analysis. We thank two anonymous 
referees for their constructive criticism and comments on our 
manuscript. 

References 

Aird D, Ross MG, Chen WS, Danielsson M, Fennell T, Russ C, Gnirke A. 

2011. Analyzing and minimizing PGR amplification bias in Illumina 

sequencing libraries. Genome Biology 12: R18. 
Alabadi D, Oyama T, Yanovsky MJ, Harmon FG, Mas P, Kay SA. 2001. 

Reciprocal regulation between TOC1 and LHY/CCA1 within the Arabidopsis 

circadian clock. Science 293: 880—883- 
Anwer MU, Davis SJ. 2013. An overview of natural variation studies in the 

Arabidopsis thaliana circadian clock. Seminars in Cell & Developmental Biology. 

doi: 10.10l6/j.semcdb.2013.03.006. 
Badr A, Mueller K, Schaefer-Pregl R, El-Rabey H, Efgen S, Ibrahim HH, Pozzi 

C, Rohde W, Salamini F. 2000. On the origin and domestication history of 

barley {Hordeum vulgare). Molecular Biology and Evolution 17: 499—510. 
Beales J, Turner A, Griffiths S, Snape JW, Laurie DA. 2007. A pseudo-response 

regulator is misexpressed in the photoperiod insensitive Ppd-Dla mutant of 

wheat (Triticum aestivum L.). Theoretical and Applied Genetics 115: 721—733. 
Borner A, Buck-Sorlin GH, Hayes PM, Malyshev S, Korzun V. 2002. Molecular 

mapping of major genes and quantitative trait loci determining flowering time 

in response to photoperiod in barley. Plant Breeding 121: 129—132. 
Campoli C, Drosse B, Searle I, Coupland G, von Korff M. 2012a. Functional 

characterisation of HvCOl, the barley {Hordeum vulgare) flowering time 

ortholog of CONSTANS. Plant Journal®): 868-880. 

©2013 The Authors 

New Phytologist© 2013 New Phytologist Trust 



Campoli C, Shtaya M, Davis SJ, von KorffM. 2012b. Expression conservation 
within the circadian clock of a monocot: natural variation at barley Ppd-Hl 
affects circadian expression of flowering time genes, but not clock orthologs. 
BMC Plant Biology 12: 97. 

Carre IA. 2002. ELF3: a circadian safeguard to buffer effects of light. Trends in 
Plant Science!: 4—6. 

Close TJ, Bhat PR, Lonardi S, Wu Y, Rostoks N, Ramsay L, Druka A, Stein N, 

Svensson JT, Wanamaker S etaL 2009. Development and implementation of 

high-throughput SNP genotyping in barley. BMC Genomics 10: 582. 
Clotault J, Thuillet AC, Buiron M, De Mita S, Couderc M, Haussmann BIG, 

Mariac C, Vigouroux Y. 2012. Evolutionary history of pearl millet (Pennisetum 

glaucum (L.) R. Br.) and selection on flowering genes since its domestication. 

Molecular Biology and Evolution 29: 1 199-1212. 
Comadran J, Kilian B, Russell J, Ramsay L, Stein N, Ganal M, Shaw P, Bayer 

M, Thomas W, Marshall D etaL 2012. Natural variation in a homolog of 

Antirrhinum CENTRORADIALIS contributed to spring growth habit and 

environmental adaptation in cultivated barley. Nature Genetics 44: 1388—1392. 
Crooks GE, Hon G, Chandonia JM, Brenner SE. 2004. WebLogo: a sequence 

logo generator. Genome Research 1 4: 1188—1190. 
Dai S, Wei X, Pei L, Thompson RL, Liu Y, Heard JE, RuffTG, Beachy RN. 

2011. BROTHER OF LUX ARRHYTHMO is a component of the 

Arabidopsis circadian clock. Plant CellTb: 961—972. 
Davis SJ. 2002. Photoperiodism: the coincidental perception of the season. 

Current Biology 12: R841-R843. 
DePristo MA, Banks E, Poplin R, Garimella KV, Maquire JR, Hartl C, 

Philippakis AA, del Angel G, Rivas MA, Hanna M et at 20 1 1 . A framework 

for variation discovery and genotyping using next-generation DNA sequencing 

data. Nature Genetics 43: 491-498. 
Dixon LE, Knox K, Kozma-Bogmar L, Southern MM, Pokhilko A, Millar AJ. 

2011. Temporal repression of core circadian genes is mediated through EARLY 

FLOWERING 3 in Arabidopsis. Current Biology -21: 120-125. 
Druka A, Franckowiak J, Lundqvist U, Bonar N, Alexander J, Houston K, Radovic 

S, Shahinnia F, Vendramin V, Morgante M etaL 2011. Genetic dissection of 

barley morphology and development. Plant Physiology 155: 617—627. 
Ehrenreich IM, Hanzawa Y, Chou L, Roe JL, Kover PX, Purugganan MD. 

2009. Candidate gene association mapping of Arabidopsis flowering time. 

Genetics 183: 325-335. 
Faure S, Turner AS, Gruszka D, Christodoulou V, Davis SJ, von Korff M, 

Laurie DA. 2012. Mutation at the circadian clock gene EARLY MATURITY 8 

adapts domesticated barley (Hordeum vulgare) to short growing seasons. 

Proceedings of the National Academy of Sciences, USA 109: 8328—8333- 
Favret EA, Ryan GS. 1966. New useful mutants in plant breeding. In: Mutations 

in Plant Breeding. Vienna, Austria: International Atomic Energy Agency, 49-6 1 . 
Gallagher LW, Franckowiak JD. 1997. Description of Stock number BGS130. 

Barley Genetics Newsletter 26: 166. 
Gallagher LW, Soliman KM, Vivar H. 1991. Interactions among loci conferring 

photoperiod insensitivity for heading time in spring barley. Crop Science 31: 

256-261. 

Hall TA. 1999. BioEdit: a user-friendly biological sequence alignment editor and 
analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series 41: 
95-98. 

Hazen SP, SchultzTF, Pruneda-Paz JL, Borevitz JO, Ecker JR, Kay SA. 2005. 

LUX ARRHYTHMO encodes a Myb domain protein essential for circadian 

rhythms. Proceedings of the National Academy of Sciences, USA 102: 10387—10392. 
Heifer A, Nusinow DA, Chow BY, Gehrke AR, Bulyk ML, Kay SA. 2011. LUX 

ARRHYTHMO encodes a nighttime repressor of circadian gene expression in 

the Arabidopsis core clock. Current Biology 21: 126-133- 
Herrero E, Kolmos E, Bujdoso N, Yuan Y, Wang M, Berns MC, Uhlworm H, 

Coupland G, Saini R, Jaskolski M etaL 2012. EARLY FLOWERING4 

recruitment of EARLY FLOWERING3 in the nucleus sustains the Arabidopsis 

circadian clock. Plant OZ/24: 428-443. 
Higgins JA, Bailey PC, Laurie DA. 2010. Comparative genomics of flowering 

time pathways using Brachypodium distachyon as a model for the temperate 

grasses. PLoS ONE5: el0065. 
Huang W, Perez-Garci'a P, Pokhilko A, Millar AJ, Antoshechkin I, Riechmann 

JL, Mas P. 2012. Mapping the core of the Arabidopsis circadian clock defines 

the network structure of the oscillator. Science 336: 75—79. 

New Phytologist (2013) 199: 1045-1059 
www.newphytologist.com 



058 Research 



New 

Phytologist 



Hiibner S, Hbffken M, Oren E, Haseneyer G, Stein N, Graner A, Schmid K, 

Fridman E. 2009. Strong correlation of the population structure of wild barley 

{Hordeum spontaneum) across Israel with temperature and precipitation 

variation. Molecular Ecology 18: 1523—1536. 
Huson DH, Bryant D. 2006. Application of phylogenetic networks in 

evolutionary studies. Molecular Biology and Evolution 23: 254—267- 
Huson DH, Scornavacca C. 2012. Dendroscope 3: an interactive tool for rooted 

phylogenetic trees and networks. Systematic Biology 61: 1061—1067- 
International Barley Sequencing Consortium (IBSC). 2012. A physical, 

genetic and functional sequence assembly of the barley genome. Nature 491: 

711-716. 

Izawa T. 2007. Adaptation of flowering- time by natural and artificial selection in 
Arabidopsis and rice. Journal of Experimental Botany 58: 309 1—3097. 

Jones H, Leigh FJ, Mackay I, Bower MA, Smith LMJ, Charles MP, Jones G, 
Jones MK, Brown TA, Powell W. 2008. Population-based resequencing 
reveals that the flowering time adaptation of cultivated barley originated east of 
the Fertile Crescent. Molecular Biology and Evolution 25: 221 1—2219- 

Katoh K, Toh H. 2008. Recent developments in the MAFFT multiple sequence 
alignment program. Briefings in Bioinformatics*): 286-298. 

Kolmos E, Herrero E, Bujdoso N, Millar AJ, Toth R, Gyula P, Nagy F, Davis 
SJ. 201 1. A reduced-function allele reveals that EARLY FLOWERING3 
repressive action on the circadian clock is modulated by phytochrome signals in 
Arabidopsis. Plant Cell23: 3230-3246. 

Kolmos E, NowakM, Werner M, Fischer K, Schwarz G, Mathews S, Schoof H, 
Nagy F, Bujnicki JM, Davis SJ. 2009. Integrating ELF4 into the circadian 
system through combined structural and functional studies. HFSP Journal^: 
350-366. 

von Korff M, Leon J, Pillen K. 2010. Detection of epistatic interactions between 
exotic alleles introgressed from wild barley (H. vulgaressp. spontaneum) . TAG. 
Theoretical and Applied Genetics. 121: 1455-1464. 

von Korff M, Wang H, Leon J, Pillen K. 2006. AB-QTL analysis in spring 
barley: II. Detection of favourable exotic alleles for agronomic traits 
introgressed from wild barley (H. vulgaressp. spontaneum). Theoretical and 
Applied Genetics 112: 1221-1231. 

Laurie DA, Pratchett N, Bezant JH, Snape JW. 1995. RFLP mapping of five 
major genes and eight quantitative trait loci controlling flowering time in a 
winter x spring barley Hordeum vidgareh. cross. Genome 38: 575—585- 

Li H, Durbin R. 2009a. Fast and accurate short read alignment with Burrows- 
Wheeler transform. Bio informatics 25: 1754—1760. 

Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, 
Abecasis G, Durbin R. 2009b. The Sequence Alignment/Map format and 
SAM tools. Bioinformatics 25: 2078-2079. 

Lindqvist C, Laakkonen L, Albert VA. 2007. Polyglutamine variation in a 

flowering time protein correlates with island age in a Hawaiian plant radiation. 
BMC Evolutionary Biology 1 : 105. 

Lou P, Wu J, Cheng F, Cressman LG, Wang X, McClung CR. 2012. 

Preferential retention of circadian clock genes during diploidization following 
whole genome triplication in Brassica rapa. Plant CelllA: lAl^—lAlG. 

Lundqvist U. 2009. Eighty years of Scandinavian barley mutation genetics and 
breeding. In: Shu QY, ed. Induced Plant Mutations in the Genomics Era. Rome, 
Italy: FAO, 39-43- 

Lyons E, Freeling M. 2008. How to usefully compare homologous plant genes 
and chromosomes as DNA sequences. Plant Journal^: 661—673. 

Mayer KF, Martis M, Hedley PE, Simkova H, Liu H, Morris JA, Steuernagel 
B, Taudien S, Roessner S, Gundlach H etal. 2011. Unlocking the barley 
genome by chromosomal and comparative genomics. Plant Cell 23: 1249— 
1263. 

Mayer KF, Taudien S, Martis M, Simkova H, Suchankova P, Gundlach H, 
Wicker T, Petzold A, Felder M, Steuernagel B etal. 2009. Gene content and 
virtual gene order of barley chromosome 1H. Plant Physiology 151: 496—505- 

McWatters HG, BastowRM, Hall A, Millar AJ. 2000. The ELF3 zeitnehmer 
regulates light signalling to the circadian clock. Nature 408: 716—720. 

Murakami M, Tago Y, Yamashino T, Mizuno T. 2007. Comparative overviews 
of clock- associated genes of Arabidopsis thaliana and Oryza sativa. Plant and 
Cell Physiology 48: 110-121. 

Murphy RL, Klein RR, Morishige DT, Brady JA, Rooney WL, Miller FR, 
Dugas DV, Klein PE, Mullet JE. 20 1 1 . Coincident light and clock regulation 

New Phytologist (2013) 199: 1045-1059 
www.riewphytologist.com 



of pseudoresponse regulator protein 37 '{PRR37) controls photoperiodic flowering 
in sorghum. Proceedings of the National Academy of Sciences, USA108: 16469— 
16474. 

Nagel DH, Kay SA. 2012. Complexity in the wiring and regulation of plant 

circadian networks. Current Biology 22: R648— R657- 
Niwa Y, Takafumi Y, Takeshi M. 2009. The circadian clock regulates the 

photoperiodic response of hypocotyl elongation through a coincidence 

mechanism in Arabidopsis thaliana. Plant and Cell Physiology 50: 838—854. 
Nusinow DA, Heifer A, Hamilton EE, King JJ, Imaizumi T, Schultz TF, Farre 

EM, Kay SA. 2011. The ELF4-ELF3-LUX complex links the circadian clock 

to diurnal control of hypocotyl growth. Nature 475: 398—402. 
Onai K, Ishiura M. 2005. PHYTOCLOCK 1 encoding a novel GARP protein 

essential for the Arabidopsis circadian clock. Genes to Cells 10: 963—972. 
Pokhilko A, Fernandez AP, Edwards KD, Southern MM, Halliday KJ, Millar 

AJ. 2012. The clock gene circuit in Arabidopsis includes a repressilator with 

additional feedback loops. Molecular Systems Biology 8: 574. 
Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R. 

2005. InterProScan: protein domains identifier. Nucleic Acids Research 'id: 

W116-W120. 

Saito H, Ogiso-Tanaka E, Okumoto Y, Yoshitake Y, Izumi H, Yokoo T, 
Matsubara K, Hori K, Yano M, Inoue H etal. 2012. ^/'encodes an ELF3- 
like protein and promotes rice flowering by negatively regulating the floral 
repressor gene Ghd7 under both short- and long-day conditions. Plant and Cell 
Physiology 53: 717-728. 

Salome PA, To JP, Kieber JJ, McClung CR. 2006. Arabidopsis response 

regulators ARR3 and ARR4 play cytokinin-independent roles in the control of 
circadian period. Plant Cell 18: 55—69. 

SAS Institute. 2009. The SAS system for Windows, release 9. 1.3. Cary, NC, USA: 
SAS Institute Inc. 

Schoning JC, Staiger D. 2005. At the pulse of time: protein interactions 
determine the pace of circadian clocks. FEBS Letters 579: 3246—3252. 

Sonnhammer ELL, Koonin EV. 2002. Orthology, paralogy and proposed 
classification for paralog subtypes. Trends in Genetics 18: 619—620. 

Stamatakis A, Hoover P, Rougemont J. 2008. A rapid bootstrap algorithm for 
the RAxML Web servers. Systematic Biology 1 )!: 758-771. 

Takata N, Saito S, Saito CT, Uemura M. 2010. Phylogenetic footprint of the 
plant clock system in angiosperms: evolutionary processes of pseudo-response 
regulators. BMC Evolutionary Biology 10: 126. 

Turner A, Beales J, Faure S, Dunford RP, Laurie DA. 2005. The pseudo- 
response regulator Ppd-Hl provides adaptation to photoperiod in barley. 
Science^: 1031-1034. 

Undurraga SF, Press MO, Legendre M, Bujdoso N, Bale J, Wang H, Davis SJ, 
Verstrepen KJ, Queitsch C. 2012. Background-dependent effects of 
polyglutamine variation in the Arabidopsis thaliana gene ELF3. Proceedings of 
the National Academy of Sciences, USA 109: 19363-19367- 

Waddington SR, Cartwright PM, Wall PC. 1 983. A quantitative scale of spike 
initial and pistil development in barley and wheat. Annals of Botany 51: 119—130. 

Wang G, Schmalenbach I, von Korff M, Leon J, Kilian B, Rode J, Pillen K. 
2010. Association of barley photoperiod and vernalization genes with QTLs for 
flowering time and agronomic traits in a DH-population and a set of wild 
barley introgression lines. Theoretical and Applied Genetics 120: 1559-1574. 

Weller JL, Liew LC, Hecht VF, Rajandran V, Laurie RE, Ridge S, Wenden B, 
Vander Schoor JK, Jaminon O, Blassiau C etal. 2012. A conserved molecular 
basis for photoperiod adaptation in two temperate legumes. Proceedings of the 
National Academy of Sciences, USA 109: 21158-21163. 

Zagotta MT, Hicks KA, Jacobs CI, Young JC, Hangarter RP, Meeks-Wagner 
DR. 1996. The Arabidopsis ELF3 gene regulates vegetative 
photomorphogenesis and the photoperiodic induction of flowering. Plant 
JoumallO: 691-702. 

Zakhrabekova S, Gough SP, Braumann I, Miiller AH, Lundqvist J, Ahmann K, 
Dockter C, Matyszczak I, Kurowska M, Druka A etal. 2012. Induced 
mutations in circadian clock regulator Mat-a facilitated short-season adaptation 
and range extension in cultivated barley. Proceedings of the National Academy of 
Sciences, USA 109: 4326-4331. 

Zhao J, Huang X, Ouyang X, Chen W, Du A, Zhu L, Wang S, Deng XW, Li S. 
2012. OsELF3-l, an ortholog of Arabidopsis EARLY FLOWERING 3, regulates 
rice circadian rhythm and photoperiodic flowering. PLoS ONE7: e43705- 

©2013 The Authors 
New Phytologist © 20 1 3 New Phytologist Trust 



New 

Phytologist 




Supporting Information 

Additional supporting information may be found in the online 
version of this article. 



Table S5 Results of the gene ontology analysis of Brachypodium 
genes surrounding the ZKflocus 



Methods SI Preparation of TruSeq libraries. 



Fig. SI Expression patterns of clock output genes. 



Methods S2 Read processing workflow using the Galaxy server 
(http : / /galaxy . wur . nl) . 



Fig. S2 Location of genes and SNPs distinguishing between 
Bowman and Bowman(eaml 0) on barley linkage groups. 



Methods S3 Polymorphism filtering workflow. 



Fig. S3 Multiple alignment of LUX-like protein sequences from 
19 plant species. 



Methods S4 Permanent web links to the detailed results of the 
synteny analysis in the vicinity of LUX genes in three Poaceae 
species. 



Table SI AN OVA for meristem development and gene expres- 
sion differences 



Notes SI Sequencing of the barley ARR-Yike gene. 



Table S2 Accessions used for the resequencing of HvLUXl 



Please note: Wiley-Blackwell are not responsible for the content 
or functionality of any supporting information supplied by the 
authors. Any queries (other than missing material) should be 
directed to the New Phytologist Central Office. 



Table S3 Barley flowering-related genes selected for targeted 
enrichment 



Table S4 PCR primers specific for the barley ARR-like gene 




About New Phytologist 



• New Phytologist is an electronic (online-only) journal owned by the New Phytologist Trust, a not-for-profit organization dedicated 
to the promotion of plant science, facilitating projects from symposia to free access for our Tansley reviews. 

• Regular papers, Letters, Research reviews, Rapid reports and both Modelling/Theory and Methods papers are encouraged. 

We are committed to rapid processing, from online submission through to publication 'as ready' via Early View- our average time 
to decision is <25 days. There are no page or colour charges and a PDF version will be provided for each article. 

• The journal is available online at Wiley Online Library. Visit www.newphytologist.com to search the articles and register for table 
of contents email alerts. 

• If you have any questions, do get in touch with Central Office (np-centraloffice@lancaster.ac.uk) or, if it is more convenient, 
our USA Office (np-usaoffice@ornl.gov) 

• For submission instructions, subscription and all the latest information visit www.newphytologist.com 



© 2013 The Authors New Phytologist (2013) 199: 1045-1059 

New Phytologist© 2013 New Phytologist Trust WWW.newphytologist.com 



